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Stability  of  a  Thick  Elastic  Plate  Under  Thrust 

by 

K.N.  Sawyers  and  R.S.  Rivlin 
Lehigh  University,  Bethlehem,  Pa. 


ABSTRACT 


When  a  rectangular  plate  of  incompressible  neo-Hookean 
elastic  material  is  subjected  to  a  thrust,  bifurcations  of  the 
flexural  or  barreling  types  become  possible  at  certain  critical 
values  of  the  compression  ratio.  The  states  of  pure  homogeneous 
deformation  corresponding  to  these  critical  compression  ratios 
are  states  of  neutral  equilibrium.  Their  stability  is  investi¬ 
gated  on  the  basis  of  an  energy  criterion,  without  restriction 
on  the  thickness  of  the  plate. 

The  critical  state  corresponding  to  the  lowest  order  flexural 

mode  is  found  to  be  stable  (unstable)  if  the  aspect  ratio 

(thickness/length)  is  less  (greater)  than  about  0.2.  Agreement 

with  the  classical  Euler  theory  is  established  in  the  limiting 
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Introduction 


We  consider  a  rectangular  plate  of  incompressible  isotropic 
elastic  material,  with  neo-Hookean  strain-energy  function,  to 
be  situated  with  its  edges  parallel  to  the  axes  of  a  rectangular 
cartesian  coordinate  system  x  .  The  plate  is  acted  on  by  forces, 
applied  as  dead-loads  normally  to  the  faces  of  the  plate  which 
are  perpendicular  to  the  1  and  3  axes  of  the  reference  system. 

The  faces  perpendicular  to  the  2 -axis  are  force-free.  The  con¬ 
straints  on  the  1  and  3 -faces  are  such  as  to  permit  the  plate 
to  undergo  pure  homogeneous  deformations  under  the  action  of 
the  applied  forces.  The  load  in  the  1-direction  is  assumed  to 
be  a  thrust  and  at  certain  critical  values  of  this  thrust 
bifurcations  occur. 

If  the  load  in  the  3-direction  is  a  tension  and  the  3- 
dimension  of  the  plate  is  sufficiently  large,  we  may  assume  that 
these  bifurcations  are  plane  strains  in  the  12-plane,  possibly 
superposed  on  a  uniform  extension  in  the  3-direction.  These  plane 
strains  may  be  of  the  flexural  or  barreling  types.  The  latter  have 
little  practical  interest  (see  §8).  They  are  discussed  in  the 
present  paper  in  the  interest  of  completeness,  since  their  in¬ 
clusion  involves  negligible  complication  of  the  discussion. 

The  critical  compression  ratios  in  the  1-direction  at  which 
bifurcations  occur  have  been  previously  discussed  by  a  number  of 
authors  with  various  degrees  of  generality.  For  flexural  bifurca¬ 
tions  of  a  neo-Hookean  incompressible  material  they  were  first 
determined  by  Biot  [1] .  They  were  determined  by  Sawyers  and 
Rivlin  [2] ,  for  an  arbitrary  incompressible  isotropic  elastic 
material,  in  both  the  flexural  and  barreling  cases. 
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In  the  present  paper  we  discuss  the  stability  of  the  states 
of  pure  homogeneous  deformation  at  which  bifurcations  occur.  The 
stability  criterion  employed  is  the  energy  criterion  and  the  pro¬ 
cedure  used  is  essentially  that  of  Koiter  [3,4,5]*.  An  equilibrium 
state  of  the  system  is  a  state  of  deformation  at  which  the  total 
potential  energy  of  the  system  has  a  stationary  value,  i.e.  at  which 
its  first  variation  is  zero.  This  state  will  be  stable  or  unstable 
accordingly  as  this  stationary  value  is  or  is  not  a  proper  minimum, 
i.e.  accordingly  as  the  second  variation  of  the  potential  energy  is 
or  is  not  positive  definite. 

A  state  of  pure  homogeneous  strain  corresponding  to  a  bifur¬ 
cation  is  one  of  neutral  equilibrium.  It  will  be  stable  if  the 
potential  energy  is  greater  for  every  state  lying  in  some  neighbor¬ 
hood  of  it,  which  satisfies  the  kinematic  constraints.  It  will  be  . 
unstable  if  the  potential  energy  is  less  for  some  such  state. 

(The  neighborhood  is  limited  to  sufficiently  small  deformations  in 
the  12-plane  superposed  on  a  uniform  extension  in  the  3-direction.) 
We  determine  whether  this  is  the  case  by  calculating  the  stationary 
value  of  this  excess  potential  energy.  The  critical  pure  homogene¬ 
ous  deformation  is  stable  (unstable)  if  this  stationary  value  is 
positive  (negative) . 

These  calculations,  which  are  rather  cumbersome,  are  carried 
out  without  restriction  on  the  magnitude  of  the  aspect  ratio 
(2-dimension/l-dimension)  of  the  plate.  They  lead  to  the  conclusion 
that,  for  the  lowest  mode  of  flexural  bifurcation,  the  pure  homo¬ 
geneous  state  is  stable,  provided  that  the  aspect  ratio  is  less  than 
about  0.2  and  is  unstable  for  higher  aspect  ratios.  The  stability 

*  In  formulating  the  problem  we  have  had  the  benefit  of  extensive  discussion 
and  correspondence  with  Professor  Koiter,  for  which  we  are  extremely  grateful. 

Most  of  his  suggestions  have  been  incorporated  in  this  paper,  which  owes  much 
to  his  generous  help  and  advice. 
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at  low  aspect  ratios  is  studied  by  means  of  an  asymptotic 
formula  for  the  stationary  value  of  the  excess  potential  energy 
which  is  valid  up  to  degree  5  in  the  aspect  ratio.  This  agrees 
well  with  the  exact  calculations  up  to  fairly  high  aspect  ratios. 
Moreover,  in  the  limiting  case  as  the  aspect  ratio  tends  to  zero, 
it  agrees  precisely  with  the  result  derived  by  Euler  on  the  basis 
of  elastica  theory. 
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2.  Statement  of  the  problem 

We  consider  a  rectangular  plate  of  incompressible  neo-Hookean 
elastic  material,  which  has  its  edges  parallel  to  the  axes  of  a 
rectangular  cartesian  coordinate  system  x  .  Let  C  be  the  vector 
position,  relative  to  the  origin  of  the  system  x  ,  of  a  generic 
particle  of  the  plate  in  its  undeformed  state  (state  0)  and  let  its 
bounding  surfaces  in  this  state  be 


51  *  ±ll  »  52  ’  ±l2  »  *3  *  ±l3  U3»l19l2)  .  (2.1) 

We  suppose  that  the  plate  is  maintained  in  an  equilibrium 
state  of  pure  homogeneous  deformation  (state  I) ,  with  extension 
ratios  X^Xg^  an<*  PrinciPal  directions  parallel  to  the  coordi¬ 
nate  axes,  by  normal  forces  applied  to  the  surfaces  ^  »  tl±  and 
?3  *  ±f3  ,  the  surfaces  £2  »  ±l2  remaining  force-free.  Let  the 
resultant  loads  applied  to  the  surfaces  ■  tt^  and  53  *  ±f3 
be  ±RX  and  ±R3  respectively  and  suppose  that  Rx  £  0  (thrust) . 

In  state  I,  we  suppose  that  the  surfaces  initially  at 
51  *  ±l1  and  C3  *  ±^3  are  constrained  so  that,  in  the  deforma¬ 
tion,  they  move  parallel  to  the  1  and  3-axes  respectively,  but 
points  on  them  are  free  to  move  normal  to  these  directions  (i.e. 
the  tangential  components  of  the  surface  tractions  are  zero) . 

Let  X  be  the  vector  position  in  state  I  of  the  particle 
which  has  vector  position  5  in  state  0  .  Then, 


Ml  ’ 


M2  ■ 


Ms  ' 


X _ X _ 
12  3 


1.  (2.2) 
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We  now  consider  that  the  plate  undergoes  a  further  deforma¬ 
tion  which  is  the  superposition  of  a  uniform  stretch  in  the  3- 
direction  and  a  plane  deformation  parallel  to  the  12-plane.  We 
call  the  state  which  is  then  reached  state  II.  We  suppose  that 
the  particle  which  is  at  X  in  state  I  moves  to  x  in  state  II, 
where 


x  ■  X  +  u  , 


(2.3) 


and 


U1  ■  ui(C;l>C2)>  u2  *  U2(51*525,  u3 


X3E53 


(2.4) 


Since  the  material  is  incompressible  u  must  Satisfy  the  relation 


(1+E) (XjUg  2+X2uT  -  + 


1,1 


U1 ,1U2 ,2 


U1 , 2U2 , 1^ 


X1X2E  " 


(2.5) 


Let  Wj  and  Wjj  denote  the  strain  energies,  per  unit 
volume,  in  states  I  and  II  respectively.  Since  the  material  of 
the  plate  is  neo-Hookean,  they  are  given,  in  appropriate  units,  by 


WI  =  I(Xi+X2+X3'3)  ’ 


W 


+  X2^2  +  XlU  +  X2*l  +  X3C1+E)" 


II 


With  (2.2),  (2.3)  and  (2.4),  we  obtain  from  (2.6) 


WII  WI  Xlul,l  +  X2U2,2  +  X3E 


+  lK!l  +  *2%  +  Ul!2  +  U2!l  +  ^ 


(2.6) 


(2.7) 
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We  consider  that  state  I  is  an  equilibrium  state.  We  con¬ 
sider  also  that  in  state  II  the  resultant  loads  on  the  faces  of 
the  plate  are  the  same  as  in  state  I  and  that,  while  state  II  is 
not  necessarily  an  equilibrium  state,  in  it  the  particles  of  the 
plate  are  at  rest,  at  any  rate  instantaneously. 

The  Piola-Kirchhoff  stress  in  state  I,  denoted  H  .  is 

ai 

given  by 


X2  X2 

nil  X1  '  17  *  n33  *  X3  "  XT  * 


(2.8) 


the  remaining  components  of  n  ^  being  zero.  Then, 


R1  "  Al2^3Rll  *  R3  “  4*1^2n33  * 


(2.9) 


It  follows  that,  in  passing  from  state  I  to  state  II,  the  in¬ 
crease  in  potential  energy  of  the  loads  is 


-  R-^u^  ^  a  £  -  ZR^X^E^ 


l  l 

r  2  ,1 


2l^j  |  (nnui,i  +  n33X3E^d^ld^2  . 


(2.10) 


-l  -i 

2 


Let  Gj  and  Gjj  be  the  potential  energies  of  the  system 
(plate  and  load)  in  states  I  and  II  respectively.  From  (2.7)  and 
(2.10),  we  obtain,  with  (2.8), 
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def . 

G  -  Gjj-Gj  -  2i3 


I^U1,1  +  U2,2  +  U1 , 2  +  U2%  +  X3E  ^  |  d^ld^2  *  C2*11) 


l0  A  (\2 

\  (xf  U!,l  *  V2,2  *  *! 

-l  -l  'o  X 
2  1 


Using  (2.5)  to  substitute  for  u  in  (2.11)  we  can  rewrite 

^  9  £ 


(2.11)  as 


rl2  tl l  f  X 


-£  -l 
c2  <-1 


°  "  Ull  I  jxf  fUl,2U2,l  '  U1,1U2,25 


XgE2 


*  Wl  *  “2?2  *  “X%  *  »2%  *  ^  «1«2  •  '2-12) 


The  necessary  and  sufficient  condition  for  stability  of  the 
plate  in  state  I,  under  the  specified  loading  conditions,  is  that 
G  be  positive  definite  for  all  u  and  E  lying  in  a  neighbor¬ 
hood  of  u  *  0  ,  E  =  0  and  satisfying  the  constraint  (2.5) 
throughout  the  plate,  as  well  as  the  kinematic  constraints 


E  =  constant  everywhere. 


ux  «  constant  on  C1  =  ±^1 


(2.13) 


We  shall  investigate  the  validity  of  this  condition  when 
state  I  is  a  critical  state  for  existence  of  a  bifurcation 


solution. 
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3.  The  bifurcation  solution 

A  necessary  condition  for  stability  of  state  I  is  that  the 
second  variation  of  G  be  non-negative  for  all  sufficiently  small, 
kinematically  admissible  values  of  u^u^E  .  This  second  varia¬ 
tion  G2[u1,u2,E]  is  given,  from  (2.12),  by 


*2  ,ll  [  X, 


G2[ui,u2,E]  ■  2£3J  |  |  X“  ^ul,2U2,l“Ul,lU2,2^  +  X2E 

-l  -l  L  1 

*  *  X1e2)]  d?ld52 


(3.1) 


u1,u2,E  must  satisfy  the  kinematic  constraint  implied  by  the 
linearized  incompressibility  condition 


ul,l  +  XJ  u2,2  +  E  "  0  ’  C3*2 

which  is  obtained  from  (2.5),  and  by  the  conditions  (2.13). 

A  sufficient  condition  for  stability  of  state  I  is  that 
G2[u1,u2,E]  be  positive  definite.  Thus,  state  I  is  at  the 
stability  limit  if  G2[u1,u2,E]  has  a  zero  minimum  for 
some  non-vanishing  ulfu2,E  . 

We  shall  now  determine  the  values  of  u1,u2,E  for  which  Gg 
has  a  stationary  value,  subject  to  kinematic  constraints  (3.2) 
and  (2.13).  These  are  given  by  setting  the  first  variation 
<SG 2  of  G2  equal  to  zero.  We  take  account  of  the  constraint 
(3.2)  by  the  method  of  Lagrange  multipliers  and  denote  the  multi¬ 
plier  of  (3.2)  by  _2£3p  .  We  then  obtain,  with  the  notation 


x  *  x2/x1  , 
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/*2  ,L  1  f 

I  J  UU1,1‘XU2S2  •  fr  )6ul,l  +  CU2,2-Xu1,1-!t)5u2,2 
“42  L 


6G2  -  U3 


+  CU1,2  +  XU2,1)5U1,2  +  <U2,1  +  XU1, 2^2,1 


+  (2X*E  +  Xp  -  p)6E  \  dC1d52  «  0 


This  relation  may  be  rewritten  as 


f  Z2  A  f  P,± 

J  ^U1 , 11+U1 , 22  ”  I” )6U1  +  CU2, 

l  1 


P » p 

11+U2,22'  X^_^<SU2 


'2  1 


- C2X|E  +  X^E  -p) 6E  ^ 


rl  2 


^Ul,l‘Xu2,2'  ^5U1  +  ^U2,1+XU1,2^  6U2^  =  _£  d^2 


1  1 


r  1 


Z2=~l2 


Equation  (3.4)  yields,  with  (Z.13), 


U1.H  +  Ul,22  •  X^P’l  =  °  ’ 


U2,ll  +  U2,22  *  X2  P,2  "  °  ’ 


t2  c1 

4£x-e2(2x2+x2)E  -  |  |  pdCxd5 


-l  -l 

2  1 


Also 


Xu2,l  +  Ul,2  "  0  » 


U2,2  '  XU1,1  ‘  $1  *  °  ’ 


on  S2  =  ±Z2 


(3.3) 


^U2,2'XU1,1-  fcSu2  +  ^Ul,2+XU2,l)<Sui]r2_  ,  =  0  *  (3-4> 


(3.5) 


(3.6) 


-11- 


and,  with  (2.13), 


(3.7) 


Equations  (3. 5) -(3. 7),  or  equivalent  equations,  have  been 
previously  obtained  by  a  number  of  workers  (see,  for  example, 

[1.2] )  using  somewhat  different  procedures.  They  yield,  with 

(3.2) ,  solutions  of  the  forms 


-sinflC1 

u' 

cos 

u  = 

1 

‘  xn  * 

U2  = 

cosfi? 

J 

sin£)5 

% 

cosfl^ 

X-B’ 

p  - 

r  ~~2  > 

CTl 

II 

O 

* 

sinfiS 

X 

where 

1  =  MT/2l1 

(n=l , 2 , . . 

the  upper  (lower) 

solution  correspondinj 

In  (3.8)  U  is  a 

function  of 

S2  only 

differentiation  with  respect 

to  52  • 

6  *  U"  - 

ft2U  , 

U  , 


(3.8) 


(3.9) 


prime  denotes 


(3.10) 


and  U  satisfies  the  differential  equation 
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u(lv>  -  (x2+i)a2u”  +  x2nuu  -  o 


(3.11) 


and  the  boundary  conditions 

U"  +  X2Q2U  »  0  , 

U,M  -  (2X2  +  1) fl2U *  -  0  , 


on  C2 


(3.12) 


We  obtain  two  possible  solutions  to  (3.11)  and  (3.12). 

With  the  notation 

n  -  Ot2  ,  (3.13) 


one  of  these  solutions  may  be  written  as 

U(S2)  =  M(coshXn52-mcoshf2C2)  , 

_  _  2X2coshXn 

m  =  — - -  , 

(X  +l)coshn 

where  X  and  n  satisfy  the  secular  equation 


(3.14) 


4X3tanhn  *  (X2+l) 2tanhXn  .  (3.15) 

The  other  solution  may  be  written  as 

U(C2)  ■  M(sinhxnc2-msinhfi?2)  , 

(3.16) 

„  _  2X2sinhXn 

m  =  — - 1 —  , 

(X^+l) sinhn 
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where  A  and  n  satisfy  the  secular  equation 

4A3tanhXn  -  (X2+l)2tanhn  .  (3.17) 

We  note  that  for  any  specified  n  ,  each  of  the  equations 
(3. IS)  and  (3.17)  yields  a  unique  value  for  X  . 

The  deformations  described  by  (3.14)  are  antisymmetric  with 
respect  to  the  13-plane,  i.e.  they  are  flexural  deformations. 

Those  described  by  (3.16)  are  symmetric  with  respect  to  the  13- 
plane,  i.e.  they  are  barreling  deformations. 

The  two  secular  equations  (3.15)  and  (3.17)  can  be  rewritten 
as 

sinh(A+l) n  .  v  (A+l) {l(X+l)2f (A-l)2}  ^  (3  18) 

sinh(X-l) n  (A-l) { (A+l) 2-X (A-l) 2} 

where  v  «  1  for  flexural  deformations  (equation  (3.15))  and 
v  =  -1  for  barreling  deformations  (equation  (3.17)). 

We  note  that  each  of  the  pairs  (3.14)2,  (3.15)  and  (3.16)2, 
(3.17)  leads  to  the  result 


m2  M  Asinh2 An 
sinh2n 


(3.19) 


By  substituting  in  (3.1)  either  the  upper  or  lower  expressions 
(3.8),  we  obtain,  by  using  (3.11),  the  result 


l3ll  (t 2  r 
X2ft2 


-l. 


(3A2  +  l)n2(UU')  '  -  (UU,U  )’  +  (U'U")'^d52  . 

(3.20) 


We  now  carry  out  the  integration  and  use  (3.12)  to  obtain 
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G 


2 


0  . 


(3 


We  conclude  that  the  critical  state  for  which  Gg  has  a 
stationary  value  is  a  state  of  neutral  equilibrium. 
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4 .  Stability  in  the  critical  case  of  neutral  equilibrium 

We  have  seen  in  §2  that  the  necessary  and  sufficient  con¬ 
dition  for  the  plate  to  be  stable  in  state  I  is  that  G  ,  given 
by  (2.12),  shall  be  positive  definite  for  all  u^u^E  lying  in 
a  neighborhood  of  ^  «  ug  ■  E  *  0  and  satisfying  the  kinematic 
constraints  (2.5)  and  (2.13).  We  now  develop  this  condition  fur¬ 
ther  in  the  case  when  state  I  is  a  critical  state  of  neutral 
equilibrium  discussed  in  §3. 

Denoting  by  u2,  E  »  0  ,  p  the  values  of  u1,u2,E,p 

given  by  equations  (3.8),  we  decompose  an  arbitrary  ultu2,E  in 
the  following  manner: 

ul  *  a^i  +  vi  »  u2  ”  a^2  +  v2  »  E  "  F  »  (4.1) 

A  A 

where  v1,  v2  is  orthogonal  to  ux,  u2,  i.e. 
l2 

j  j  "  <“lV52V2)d?ld?2  '  0  ' 

-l  -l 
2  1 

and  a  is  given  by 
l2 

|  I  tS1ul«S2u2)d?1d52 

-  -i  - - 

a  *  — — - - 

f2  f41 

-I  -l 

Substituting  from  (4.1)  in  (2.12)  and  using  (3.1),  (3.5)- 
(3.7)  and  (3.21),  we  obtain 


(4.2) 


14.3) 
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M*  I1 


+  X(vl,2V2,l'Vl,lV2,2 


A2F2 

3  +  T+7- 


+  aP(X7  vl,l+  rj  v2,2>  d*ld52 


(4.4) 


where  p  is  given  by  (3.8).  From  the  incompressibility  con¬ 
dition  (2.5)  and  (3.2),  it  follows  that  v1,  v  ,  F  must  satisfy 
the  relation 


A  A 


X1X2F  2  *  ..  „ 

A  v  0  +  A  v  +  -  +  a  (u  u„  -u  u  ) 

1  d,d  d  X,1  1+p  X,X  d,d  L,d  d, 1 


+  a^u!  1V2  2  +  u2  2V1  1  "  “l  2v2  1-  ‘  ^2  1V1  2-^ 


+  v  V  -V  V  =0 

1,1  2,2  1,2  2,1 


(4.5) 


From  (3.8)  we  obtain,  by  carrying  out  the  integration  with 
respect  to  5  , 


r  2  t  1 

j  P  ^1 , 1U2 , 2  ~U1 , 2^2 , 1^  d^ld^2  *  0  * 

_  (4.6) 


-i  -l 
*~2  1 


Then,  using  (4.5)  to  eliminate  the  term  in  (4.4)  which  is  linear 
in  v  and  v  ,  we  obtain  with  (4.6) 
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G  -  u3  \  2  \  1(7(Vi!i+V2!2+Vi!2+V2!i+X3F2) 
-l2  -lx  l 


*22f2 


+  X(V1,2V2,1'V1,1V2,2)  +  — 


2  ■»  A  /  A  A  A  A  . 

a  X  p(u  -V-  9+u  v  -u.  0v  — u0  -v.  .) 

j  * j  ^  ^  2  j  &  1 j  X  X  y  2  2  9  X  2  9  X  1)2 


-  ^3P(v1#1v2i2-vlt2v2il)  dqd52 


It  follows  from  (2.13),  (4.1)  and  (3.7)  that 


(4.7) 


vl,2  *  0  £or  *  *^1  • 


(4.8) 


For  a  fixed  small  value  of  a  ,  we  now  determine  the  values 
of  v  ,  v2,  F  ,  satisfying  the  kinematic  constraints  (4.S)  and 
(4.8)  and  the  orthogonality  condition  (4.2),  for  which  G  has  a 
stationary  value.  Since  these  values  must  be  0(a2)  ,  we  intro¬ 
duce  the  notation 


v1  *  a2u1  ,  v2  ■  a2u2  ,  F  *  a2E 


(4.9) 


Then,  neglecting  terms  of  higher  degree  than  the  fourth  in  a  , 
we  obtain  from  (4.7), 

G  -  a'-G  -  ak2*3  {  *  f  *  f 

-l  -l  l 
C2  ^1  ^ 

_  _  _  _  o_o 

*  AfUl,2U2,rUl,lU2,2)  *  A2E 

2*^2, 2“l,l"’‘l,2“2,l"“2,l“l,  2^1  d5ld^2  '  t4'10) 
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Neglecting  terms  of  higher  degree  than  the  second  in  a 2  ,  the 
incompressibility  condition  (4.5)  can  be  written  as 


X  U2,2+  E  +  X3^1, 1^2, 2'^1, 2^2,1^ 


(4.11) 


and  the  orthogonality  constraint  (4.2)  may  be  written  as 


a 


2 


(u1u1  +  u2u2)dC1dC2  *  0  . 


(4.12) 


We  determine  the  values  of  u  ,  u2,  E  for  which  G  has  a 
stationary  value  by  equating  its  first  variation  to  zero,  re¬ 
laxing  the  constraints  (4.11)  and  (4.12)  by  the  method  of  Lagrange 
multipliers.  We  denote  the  Lagrange  multipliers  for  the  con¬ 
straints  (4.11)  and  (4.12)  by  -2£3a2p  and  Zl^a2x  respectively, 
where  p  is  a  function  of  2  and  x  is  a  constant.  We 

thus  obtain  from  (4.10),  after  some  algebraic  manipulation, 


f  f*2  A  _ 

H  J  <VU1  +  VU2  '  <t<SE)d?id?2 

l  Un  -l. 


<5G  =  21. 


"2  '"l 


-l  L 


+  Og^Ug 


I  l 

I  *  t2«s2] 


1 


d^2  "  0  ’ 


(4.13) 


where  6G  ,  6ux,  <5u2,  5E  denote  the  variations  of  G,  ux,  ug,  E 


and 
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« 

0, 


Ul,ll  *  Ul,22 

'  3”  p»l  “  X3^P»1^2,2 

A  A  _ 

-  XUX  , 

U2,ll  *  “2,22 

•  ^P.a  ’  X3^P»2"l,l 

A  A 

'  P»1U1,2) 

-  XU2  , 

-p  *  (A*  ♦  2a| 

)H  , 

Sl,2  *  >S2,1  * 

A  A 

X3PU2,1  * 

*2,2  -  A*l,l  - 

1  _  A  A 

XJ  p  ■  A3pul,l  * 

“l,l  -  A*2,2  - 

1  A  A 

j-  p  -  X3pu2>2  , 

“2,1  *  A“l,2  * 

A  A 

X3pUl,2  * 

(4.14) 


From  (4.13)  we  obtain,  bearing  in  mind  that  5E  is  constant 
throughout  the  plate, 

Ln  l. 


$1  =  0  ,  $2  =  0 


f  2  f  1 

,  j  j  *d51d52  -  0  .  (4.15) 


-l  -l 
2  1 


We  also  obtain,  with  (4.8),  (4.9)  and  (3.7) 

01  *  02  *  0  for  *  ±l2  , 

and 


U2  i  -  0  ,  ux  ■  ,  say  , 

r*2 

J  *1*2  *  °  ’ 


(4.16) 


for  -  ±£1  ,  (4.17) 


where  e  is  a  constant. 


Equations  (4.15) - (4.17) ,  together  with  the  conditions  (4.11)  and 
(4.12),  yield  the  values  of  uL,  u2,  E  ,  subject  to  the  kine¬ 
matic  and  orthogonality  constraints,  for  which  G  has  a  station¬ 
ary  value.  We  discuss  the  solution  of  these  equations  in  the  next 


With  the  boundary  conditions  (5.4)  and  (4.17),  we  obtain 
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the  following  solution  of  equations  (5.1)^  2  an<*  (5.2),  for 
which  u  t  u2  is  orthogonal  to  u^  ug  : 

-  Ux  sin2nci  +  , 

u .  «  U  cos2fi£;  +  V  , 

_  (5.5) 

p  *  P  COsZfl^  +  Q  , 

X  -  0  , 


where  V  and  Q  are  given  by 

V  =  1  uu'  -  *a(e  ♦  E)52  , 

2  _  (5.6) 

Q  =  i(UU')  *  +  -4-3-  B’U’  -  X2(2e  +  E) 

1  2X2£T  2 

and  Ulf  U  and  P  are  functions  of  ?2  only  ,  which  satisfy 
the  ordinary  differential  equations 


-  4fl2U1 

+  , 
X1 

0) 

zxx2& 

(s"u  - 

eV)  , 

u” 

-  4Q2U  - 

P’ 

A2 

u> 

2x2xgn2 

(e’V 

-  s’u”)  , 

(5.7) 

tF 

*  ZXQUy 

oja 

. . 

*  ■  - 

■  . - . —  _ 

and  the  boundary  conditions 


on  52  *  ±Z2  .  (5.8) 

bV  , 


U1 
U  ’ 


2Xnu  s  uxpr  B  u  • 


2xnu  .  p 

X2 


0) 


2x2x2n2 


From  (5.7) 


we  obtain 


U"  +  4X2fi2U  =  -  yj-  n2  (7X2  +  1)UU '  , 

on  C2  =  ±£2  .  (5.12) 

tT"-  4  (2X2+1)£22u'  *  —  (2X2-l)fi2u'2  , 

*2 

We  note  from  (3.14)  and  (3.16)  that 

u'2  =  y  M2n2{X2cosh2Xn?2  +  m2cosh2^52  -  v(X2+m2) 

-  2Xm[cosh(X+l)nC2  -  vcosh(X-l)ft£2] }  , 

UU*  =  \  M2n{Xsinh2XnC2  +  m2sinh2ftC2  (5.13) 

-  m[  (X+l)  sinh(X  +  l) nc2  +  v(X-l)sinh(X-l)ft£2] }  . 
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It  can  easily  be  verified,  using  (5.11),  that  equation 
(5.10)  has  a  solution  of  the  form 


U  =  M1sinh2XQC2  +  M2sinh2ft£2  -  M^sinh  (A  +  l)  Q£2 
where,  provided  that  A  ^  3  , 


M1+sinh(A-l)fiC2  , 
(5.14) 


M,  .“#l2  (illl  , 

J  4A2  (3A+1) (A  +  3) 

(5.15) 

M,  =  vmM2  tfrV  . 

4A2  (3A-1) (A-3) 

The  integration  constants  M1  and  M2  are  determined  from  the 
boundary  conditions  (5.12),  with  (5.13),  as 


2 

M  =  —  {4A2f1cosh2n  -  (A2+l)  f „sinh2n }  , 

1  32  AAA  1  * 

/  (5.16) 

=  —  {  Af  sinh2An  -  (A2+l) f_ cosh2An}  , 

16A2A  2  1 

where 

A  *  4A3sinh2 Ancosh2n  -  (A2+l) 2cosh2 Ans inhZn  , 

f  =  m{ f (A) sinh(A+l) n  -  vf (- A) sinh(A-l) n} 

-  (7A2  +  1)  (AsinhZAn  +  m2sinh2ri)  ,  ^ 

f2  ■  m{g (A)cosh(A+l) n  +  vg (-A)cosh(A-l)n} 

-  4 (2A2-1) { A2cosh2An  +  m2cosh2n  -  v(A^+m2)}  , 

with  f ( A )  and  g ( A )  defined  by 


f  (X)  =  (A+l)  {  3(Xg-l)(5XS2X*l)  +  7x2+  1}  > 

(3A+1) (A  +  3) 

g  (X)  =  5-^-x2-1^'1'-1^2  (7X2-2X+3)  +  8X(2X2-1)  . 

(3 A+l) (  A  +  3) 


(5.18) 


The  constants  E  and  e  can  be  obtained  from  the  relations 
(5.1) j  and  (4 .17) ^  which  express  the  dead-loading  conditions  on  the 
faces  and  £ 1  =  ±l1  respectively.  Using  (5.5)3  to 

substitute  for  p  in  (5.1)3,  we  obtain  with  (5.6)2, 


(3A2  +  A2)E  ♦  2A2e  -  -  \  H  , 


(5.19) 


where 


H  = 


2A2n2£ 


(u"2  +  ft2U ' 2) d£  . 


2  -l. 


(5.20) 


Substitution  from  (5.9)  and  (3.8)  in  (4.14)^  yields 


4A8‘ 


U"  +  4A2fl2U  +  JS-  (a '  +  28  *U  )  I  cos2ne_ 

ZX2  J  1 


+  2X2E  +  (3A2+ A2)  e  +  -~p  (U"2+  ^2u'2) 

Ai  _  2  21  x2n2 


u'(U"+  X2n2U) 


x,x¥ 


(5.21) 


From  (4.17)3  and  (5.21)  we  obtain,  with  (5.12),  (5.3),  (3.10) - 
(3.12)  and  (5.20), 


2 A2E  +  (3 X2  +  X  2) e  =  -  H  . 


(5.22) 


Equations  (5.19)  and  (5.22)  yield 
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_  X.H  .  _  X.XH  , 

E  =  2k  CX  _1)  ’  e  -  -  (2X+X^)  ,  (5.23) 

with 

k  -  5X3  +  3X  (XX^  +  l)  +  X^  .  (5.24) 

With  (3.14),  (3.16)  and  (3.19),  equation  (5.20)  yields 

H  =  i  M2f22(X2-l)fv  +  1  inh2Xn  \  _  (5.25) 

1  \  x*+i  2Xn  / 
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6 .  Stationary  value  of  G 

We  now  use  C S . 5 )  2  and  (3.8)  to  substitute  for  u2, 

ulf  u2,  p  in  (4.10)  and  carry  out  the  integration  with  respect 
to  £i  .  We  obtain 

l. 


G  -  11 


1^3  {  ^gl 


+  ?2  ^  2  * 


(6.1) 


where 


2 1 


X2e2  +  (2X2+X2)E2  -  e(2X  v'  +  bV)  +  v' 

1  23  2  x2tr 


*  — Vt  v'e’u'  » 

x2x2jt 


g2  =  |  {4fl2U2  +  U^2+  u'2+  4Q2U2-  4XQ  (UU1)  ' 


(6.2) 


u)3 


x2x¥ 


(uV-  2xnu'u1  +  2u"u  -  xnuu^)}  . 


By  using  (5.6)1,  (S.19)  and  (5.22)  in  (6.2)1  and  (5.9),  (5.3), 
(5.10)  in  (6 . 2) 2 ,  together  with  (3.10)  and  (3.11),  we  obtain 


2e  +  E  {.."2 .  n2IT'2_  1  \202 


{u"2+  n2u'2-  £  x2n2H  -  [u'(uVsj2u)]'} 


x2q2 


+  -ig.  (UuVtCUUV  +  bV]  , 


4X 


X^1 


(6.3) 


3w 


{ 3(u"+  X2n2U)  -  2bV}u' 


I6x2x2n2 


+  - 1 -  {a '  2+  4fi2a2  +  4a(BV-  X2n2BU) } 

32X2X2f22 

+  ,i  {u’(u’  +  4 x2n2u)  -  u[u  "  -  4n2(2X2+i)u’] 

8Xfi 


-28- 


+  “u'(3UU,m-  u'u"  -  2G2UU’) 

A2 

+  U[6u"2+  3ft2(X2+l)UU &2U ' 2 (5X2+11) 

+  3Q2(X2-1)ci  +  66 ’ U '  +  (2U"-  3B)(u"+  X2fi2U)] 


1  f  J 

*  -4  ct8  U}  . 


In  deriving  the  relation  (6.3)2*  the  following  identity  is  used: 


{n2(X2-l)a-  6  (U"+  X^U)  +  2  6  *U  * } '  =  0 


It  P  p 

T  A  <- 


»Tf  1  1  * 


(6.4) 


We  now  substitute  from  (6.3)  \in  (6.1)  and  use  the  boundary  con- 

\ 

ditions  (3.12)  and  (5.12),  together  with  (5.23),  to  obtain 


g  =  -  4£1£2£3(r1+  r2  +  r3  +  ru  -  r5)  , 


(6.5) 


where 


(7X2  +  4XX3+  1)H2  , 


\ —  [uu’(u,2+  xW)]£  , 

^  o  *2^2 


4X^2 


16  X2Xe:£2 


[-(7X2+1)UU'U'+  4 (X2+1)U ' 2U] 


(6.6) 


32X2X2S12£2  J_£ 


[ b (u"+  x2n2u)  -  2 g 'u 1 ] u 'd^2  , 


1 -  2  [a '  2  +  4n2a2+  4a(g,u'-  X2fi2gU) 

\2n2i„  o 


64X2X2fi2£2  l£ 


+  166 ' U ’ (UU ' )  ’  +  8X2n2(UU') '2]dC 


and  k  and  H  are  given  by  (5.24)  and  (5.25). 

We  now  introduce  into  (6.6)  the  expressions  (3.14),  (3.16) 
and  (5.14)  for  U  and  U  ,  together  with  (3.10),  (5.11)  and 
(5.25),  and  after  carrying  out  the  integrations,  we  obtain,  with 
the  notation 

M  «  M/2 l2  ,  (6.7) 

the  following  expressions: 


-  MUnUA  (7X2+4XX3+l)(X2-l)2[v  ♦  ^xl~]  2  ’ 

<  3  3  xd  +  l  4ATi 


M  p3X  (X  -1)  22  2 

. . .  x - sinh2Xr|{v(X  -l)m  +  2X  cosh2Xp 

2X2 (X2+l) 

+  (X2+l)m2cosh2n 


-  2Xm[ ( X+l) cosh (X+l)n+v(X-l) cosh (X-l)n] } 


MUn3A 


r,.  = 


- — I -  (8X(X2  +  1)2y1y.+  (7X2  +  1) (X2-l)2Y2sinh2Xp)  , 

128X5(X2+1)  1  J 

(6.8) 

3MUn3X-(X2-l)  3 

3  —  (ncu  +  l  uu(k,4-k)  ♦  c4(k,k-4)l 


64X 


F 


k-l 


♦  l  [s4(k-l ,3-k)  +  c4(k,k-2)]>  , 
k-l 

V  ,  +  l  Uc(k,4-k)  +  c5(k+l,k-3)l 

5  32X3  5  k-0  5  ? 

2 

+  l  [c  (k-l,3-k)  +  c  (k,k-2) ] }  , 
k-l  5 


where  the  5's  are  given  in  the  appendix  (§  9  ), 
by  (5.24)  and  yx,  Y2»  Y3  are  defined  by 


K 


is  defined 
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y  *  F1sinh2Xq  +  F2sinh2n 

-m[F3(X)sinh(X+l)n  +  vF3(-X)sinh(X-l)n]  , 

Y2  *  2XF1cosh2Xn  +  2F2cosh2n 

-m[(X+l)F3(X)cosh(X+l)n  +  v(X-l)F3(-X)cosh(X-l)n]  ,  (6.9) 

Y3  =»  X2cosh2Xn  +  m2cosh2ri  -v(X2+m2) 

-2Xm[cosh(X+l) n  -vcosh(X-l)n]  , 


with  F  ,  F2  and  F3(X)  defined  by 


F,  »  [4X2f,  cosh2n  -  (X2 

1  A  J- 

F2  “  X  tXf2sinh2Ar^  -  (X2 


F  (X)  *  24X 


+I)f2sinh2n]  , 
+l)f 1cosh2Xn]  , 


(6.10) 


Here,  A  ,  f1  and  fg  are  given  by  (5.17)  and  m  is  given  by 
(3.14)  if  the  deformation  is  flexural  and  by  (3.16)  if  it  is  of 
the  barreling  type.  In  either  case  m2  is  given  by  (3.19). 

Introducing  (6.8)  into  (6.5)  we  obtain  an  expression  for  G 
in  the  form 


G  -  4£1£^3M1*n3X3(g2-  Kgx)  , 


where 


K  .  1  (7X2+  4XX3  ♦  1)  , 


*1 


n(X2-l) 2 


X2-3 

j7+1 


sinh2Xn 
~Un - 


(6.11) 


(6.12) 


2X2g0  =  ^  ~1)--  sinh2Xn{v(X2-l)m2+2X2cosh2Xn+(X2+l)m2cosh2n 
2  X2+l 

-2Xm[ (X+l)  cosh  (X+l)rt+v(X-l)  cosh  (X-l)n] } 

- 1  -  {8X(X2^1)2y,Y^(7X2+1)  CX2-l)2Y9sinh2Xn} 

64X3(X2+1)  13  2 

-  ~^2~1  (nsu  *  I  Uu(k,4-k)  +  c,(k,k-4)] 

32X2  4  k-1  4  4 

2 

♦  l  [5u(k-l,3-k)  ♦  C  (k,k-2)]} 
k»l 

+  T^X  (n?5  U5(k,4-k)  +  c5(k+l,k-3)J 


2 

♦  l  U-(k-l,3-k)  ♦  r  Ck,k-2)]}  . 


In  §8  we  show  the  numerical  dependence  of  G  on  r\  . 
However,  before  discussing  this  we  shall  consider,  in  the  next 
section,  the  asymptotic  case  when  n  is  small,  i.e.  when  the 
aspect  ratio  sma^»  anc*  t*le  bifurcation  is  of  the 

flexural  type. 


7. 
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The  asymptotic  case  of  small  aspect  ratio 
In  principle  we  could  find  an  expression  for  G  ,  for  a 
flexural  bifurcation,  in  the  limiting  case  when  n  <<  1  »  direct¬ 
ly  from  the  expression  (6.5),  where  the  T’s  are  given  by  (6.6), 
by  expanding  the  T's  in  powers  of  n  and  neglecting  all  but  the 
leading  term  in  the  expression  for  G  so  obtained.  However,  this 
leading  term  would  be  of  degree  eleven  in  n  and  the  labor  in¬ 
volved  in  obtaining  expansions  of  this  degree  for  each  of  the  T's 
would  be  excessive.  We  recognize  that  this  difficulty  arises  from 
the  fact  that,  for  n  <<  1  ,  the  expression  (3.14)  for  U  is 

O  —  ll 

0(n  )  and  the  expression  (5.14)  for  U  is  0 (n  )  .  Accordingly, 
in  this  section  we  return  to  the  equations  (3.11)  and  (3.12)  for 
0  and  obtain  a  solution,  in  the  form  of  a  power  series  in  n  , 
which  is  normalized  so  that  u2(^1,0)  *  1  .  With  this  expression  for 

A 

U  ,  we  solve  (5.10)  and  (5.12)  to  obtain  a  corresponding  expression 
for  U  ,  in  the  form  of  a  power  series  in  n  .  With  these  ex- 

A  *y 

pressions  for  U  and  U  ,  we  then  obtain  from  (6.5)  and  (6.6) 
an  expression  for  G  ,  again  as  a  power  series  in  n  . 

We  define  the  dimensionless  thickness  coordinate  t  by 

t  -  %2/l2  .  (7.1) 


Then,  with  the  notation 


U(t)  =  U(*2t)  , 


(7.2) 


we  obtain  from  (3.11)  and  (3.12) 


7  17 


X2.l  d2U 


2  A 
+  XU 


(7.3) 


n2  dt2 


+  X^U  =  0  , 


for  t 


(7.4) 


l  d3U  ,7x2+1,i 

7  57  '  w 


The  critical  value  of  X  at  which  a  flexural  bifurcation 
can  occur  is  related  to  n  by  (3. IS).  It  has  been  shown  [  6] 
that  for  n  <<  1  ,  this  relation  may  be  written  as 


2  k  6^,8^ 

1  +  ain  +  a2n  +  a3n  +  0(n  ) 


(7.5) 


where 


ai  " 


a  -  16 

a2  TUT  ’  * ‘ * 


(7.6) 


With  (7.5),  we  can  obtain  a  solution  of  (7.3),  with  the  boundary 
conditions  (7.4),  in  the  form 


U(t)  -  l  +  A11n2t2  +  nU(A21t2  +  A22tU) 

+  n6(A31t2  ♦  A32tU  +  A33t6)  +  0(n8)  , 


(7.7) 


where 


1a  1a 

1  ’  A2i  *  7  ’  A: 


^31  "  T? 


,  A 


32*  "  17  *  A33 


(7.8) 
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We  note,  from  (3.8),  that  we  have  normalized  the  solution  (7.7), 
so  that  the  displacement  u2(£1,0)  3  ±  1  ,  the  +(-)  sign  being 
applicable  if  ^-n  ,  or  j(n-l),  is  even  (odd). 

With  (7.1)  and  the  notation 


Oct)  -  uu2t)  .  s  -  ^  0  -  (^)2}  , 


we  can  rewrite  (5.10)  and  (5.12)  as 


1  d^ 

dt 


4(x2+i: 


.  m20 

dt2 


(X2-l) 

2X2n  J  at 


\  ♦  4X2U  3  -  jgL  (7X2  +  1)U  ^  , 

n2  dt2  •  2 ^ 


1  d  U  4  f  71^  + ^ 

4t2x  .i)  t2X  l)  Urj  ■ 


With  (7.5)-(7.9),  equations  (7.10)  and  (7.11)  become 


1  d^fi  .  4(>2.1)  dfu  +  l6i25  .  8(0(Uy fn t 
n2  dt2  J 

+  n3(^-t+yt3)  +  0(n5)J 


"n^  dt 


n  dt‘ 


1  +  4x2u  =  wflX^n/4  +  in2  +  +  0(n6))  , 

n2  dt2  3  V.  3  ii5  J 

_i_  d^u  -  4(2X2+i)i  *  u>nx^n2f2  +  4n2  +  o(nu)j 

n3  dt3  ^ 


for  t  3  ±  1 


(7.9) 


(7.10) 


for  t  3  ±  1  .  (7.11) 


(7.12) 


(7.13) 


-  ** 


-35- 


It  follows  from  (7.12)  and  (7.13),  with  (7.5)  and  (7.6), 
that  U  may  be  expressed  in  the  form 


u  -  uax^tB^nt  +  n3(B21t  +  B22t3)  +  n5(B,nt  +  B„t3+  B„t5) 


3  11 


31  32  33 


+  n^B-.t  +  ...+  B..tT)  +  0(n9)  >  , 


where 


22 

1 

2  ’ 

B?1 

23 

71  ’ 

2 

BU 

1 

U3 

9  ’ 

=  TS  ‘ 

We  shall  not  require  the  values  of  B^1,  B^2  . 


(7.14) 


(7.15) 


With  (7.1),  we  now  rewrite  (6.6)  as 

ri  =  -A-  (n2+  4AA3  + 1}^2 » 

1  16X4k  3 

x~,£2  r /s  jit  r / \ 2  'l 


cox^n3 


r* [(«)** 

r3  *  7^757  -  tn2+1)5  w  at},.!  ■ 

i;  {^y*  4^s2*  ^  &  -  m*) 

♦  16  (0  ffl)*  S^fA  ffi  fflfu.  , 


x3sr  <-1 

64X3n4  ml 


.  ,,  d0  dU  d  fr,  du\  +  o  i  2  ~d  /{}  dU\~j 2 ") 

+  16  at  3t  ar  Vu  a?/  8A  |ar(,u  drj  ) 


(7.16) 


where  a  is  given  by  (7.9)  and  6  and  H  are  defined  by 
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(7.17) 


We  now  substitute  from  (7.5),  (7.7)  and  (7.14)  in  (7.16) 
and,  after  carrying  out  the  indicated  integrations,  we  obtain 


2  = 


x3« 


k  r 


nk  fl  3  _2  L  /91  ±  1  \  ^ 

[t-t”  +  5J 

4 +  Is  ^  +  0Cn6)  J » 

1+1  1  +  w  -  77  ^  + 


r,  = 


r„  * 


x3« 


r  ^  3  2  ,  1 7 

X3f*  'H>  ^  +  7?"  ^  +  0(n  ) 


x3n 


1  3  2  199  b  nr  6> 

7  ‘  4  n  +  TM  n  +  °fn  ^ 


(7.18) 


With  (7.18),  we  obtain  from  (6.5)  and  (4.10) 

U, 


C  aUC  -  a  ^l£3X3  n6 
G  ■  a  o  *  - ~ -  n 

«2 


fl*  +  0(nU) I  . 


) 


■l-K16*sb) 


(7.19) 


We  can  compare  this  result  with  a  result  which  was  derived 
by  Euler  on  the  basis  of  the  elastica  theory.  In  order  to  do  this, 
we  first  calculate  the  energy  associated  with  the  flexural  defor¬ 
mation  au1,  aUg  when  the  applied  load  is  zero.  Accordingly, 
in  (2.12)  we  take  ■  *2  -  1  and  substitute  ui»  u2’  E  = 
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au1,  au2,  0  , 
in  (3.8)  with 
value  of  G 


G  * 


where 
U  »  U 
by  G  , 


8a2l1l3 

rz — 


u  ,  u2  are  the  values  of  u^  ug  given 
given  by  (7.7).  Denoting  the  resulting 
we  obtain 


k 


n 


— 


l  +  0(n2) 


(7.20) 


Then  from  (7.19)  and  (7.20)  we  obtain,  with  *  1  , 


G 

G 


aV 
i6  d 


1  + 


o(n2) 


(7.21 


For  the  flexural  mode  of  lowest  order,  we  have  n  *  , 

and  (7.21)  may  be  rewritten  as 


(7.22 


in  agreement  with  the  result  of  Euler. 
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8.  Numerical  results 

The  critical  value  of  X  at  which  a  bifurcation  occurs  is 
related  to  n  by  (3.15)  or  (3.17),  accordingly  as  the  bifurcation 
is  of  the  flexural  or  barreling  type.  These  relations  are  plotted 
as  Curves  I  and  II  in  Fig.l.  Similar  curves  have  been  shown 
previously  by  a  number  of  authors  (see,  for  example,  [2]). 

The  approximate  relation  (cf.(7.5)  with  (7.6)) 

x  =  l  +  |  n2  +  ,  (8.1) 

valid  in  the  flexural  case  for  small  n,  is  plotted  as  Curve  III. 

The  asymptotic  expression  (7.19)  for  G  ,  valid  for  flexural 

deformations,  has  been  calculated  on  the  basis  of  an  expression 

for  u  ,  normalized  so  that  its  component  Ug(£  ,0)  *  ±  1  , 

1  1 

accordingly  as  n  ,  or  y(n-l),  is  even  or  odd.  In  order  to 
compare  with  this  expression  for  G  ,  the  exact  expression  given 
by  equations  (6.11)  and  (6.12),  we  have  to  normalize,  in  the  same 

A 

manner,  the  expression  for  u  on  which  it  is  based.  From  (3.8) 
and  (3.14)^,  we  have 

u ,,(^,0)  =  ±  M(l-m)  ,  (8.2) 

where  m  is  given  by  (3.14)2  and  the  +(-)  sign  is  taken 
accordingly  as  jn  ,  or  j(n-l)  ,  is  even  (odd).  Accordingly, 
in  order  to  normalize  our  results  in  the  desired  manner  we  must 
take 

M  “  ’Pm  and  ^  *  2Z2Tl-m)"  •  ^  (S,3) 
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With  this  expression  for  M  , 
(6.11)  as  (cf.(7.19)1) 


G 


lll3X3 

6  ll 


n  6r* 
n  g 


where 


is  given  by 


we  can  rewrite  equation 


(8.4) 


G 


* 


3 

2n3(l-m)U 


(g2‘Kgl)  , 


(8.5) 


and  m  is  given  by  (3.14) 2,  while  g  ,  g2  and  K  are  given 
by  (6.12)  with  v  =  1  . 

The  relation  between  G*  and  n  provided  by  (8.5),  with 
X^  *  1  ,  is  plotted  as  Curve  I  of  Fig. 2  for  n  <  0.6  .  Curve 
II  shows  the  relation  between  G*  and  n  provided  by  the  asymp¬ 
totic  formula  (7.19) 2,  with  *  1  ,  and  we  note  that 
the  agreement  is  good  at  the  lower  values  of  n  .  In  Fig.  3  the 
relation  between  G*  and  n  ,  provided  by  the  exact  formula 
(8.5),  with  X3  =  1  ,  is  plotted  for  a  wider  range  of  values  of  n 

Calculations  for  larger  values  of  n  than  those  covered  in  Fig. 3 

*  A 

indicate  that  the  G  vs.  n  curve  continues  smoothly,  with  G 
tending  to  -00  as  n  -*•  «  (i.e.,  X  3.383).  We  note  (see  §5) 
that  our  calculations  have  excluded  the  case  when  X  =  3  (i.e., 

n  c  1.63).  It  was  not  considered  worthwhile  to  investigate  this 

* 

case  separately,  since  our  computations  indicate  that  the  G  vs. 

n  curve  passes  smoothly  through  this  point. 

* 

From  Fig. 2  or  3,  we  see  that  G  is  positive  for  values  of 
n  below  about  0.32.  For  the  lowest  flexural  mode  this 
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corresponds,  from  (3.9)  and  (3.13)  with  n  *  1  ,  to  an  aspect 

* 

ratio  £2/^i  ~  °*20,  For  lar8er  values  of  n  ,  G  is  negative. 
Accordingly,  state  I,  the  homogeneous  state  corresponding  to  the 
critical  value  of  X  at  which  a  flexural  bifurcation  occurs,  is 
unstable  if  n  ^  0.32,  and  stable  if  n  <  0.32  and  this  stability 
is  not  preempted  by  the  appearance  of  a  mode  of  lower  order. 

We  note  that  for  a  specified  value  of  n  ,  the  critical 
value  of  X  ,  given  by  (3.15),  is  independent  of  X^  .  It 
follows  from  (3.14)2,  (5.17)  and  (5.18)  that  6,  f1  and  f2  are 
independent  of  X^  .  Then,  from  (6.9),  (6.10)  and  equations 
(9.1) -(9.4)  in  the  Appendix,  it  follows  that  y2,  y3  and  the 

S's  in  (6.12) 3  are  independent  of  X^  .  Consequently,  g1  and 
g2  are  independent  of  X3  .  Also,  for  a  fixed  value  of  X  ,  the 
quantity  K  ,  defined  by  (6.12)^,  decreases  monotonically  with 
increase  of  X3  .  Hence,  from  (8.5),  G*  changes  monotonically 
with  X3  .  However,  the  dependence  of  K  ,  and  hence  of  G*  , 
on  X3  is  very  slight. 

For  barreling  deformations,  the  exact  expression  for  G  is 
given  by  (6.11)  and  (6.12),  with  v  *  -  1.  We  cannot  normalize 
the  displacement  u  on  which  it  is  based  in  the  same  way  as  we 
did  for  flexural  deformations,  since  u2(f1,0)  *  0  for  barreling 

/v  ^ 

deformations.  Instead,  we  now  normalize  u  so  that  u 2(£1,f2) 

*  ±  1  accordingly  as  j  n,  or  y(n-l),  is  even  or  odd  by  taking 

M  =  1 _ 

2£2  ~  21 2  (sinhXn-msinhri)  ’ 


M 


(8.6) 
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where  m  is  now  given  by  the  expression  (3.16)2,  appropriate  to 
barreling  deformations.  With  this  expression  for  M  we  can  now 
rewrite  (6.11)  as 


where 


-  * 

G  »  1  |  3  G  , 

4  ll 


n3Cg2-Kg1) 
(sinhAn-msinhn) ^ 


(8.7) 


(8.8) 


G*  is  plotted  against  n  in  Fig. 4  for  the  case  when 
A^  -  1  .  As  in  the  case  of  flexural  deformations  the  dependence 
of  G*  on  A^  is  very  slight.  We  note  that  G*  ,  and  hence 
G  ,  is  positive  for  all  n  .  Accordingly,  the  homogeneous  state 
corresponding  to  a  critical  value  of  A  at  which  a  barreling 
bifurcation  occurs  is  stable. 

We  have  included  the  discussion  of  barreling  bifurcations 
for  completeness.  However,  it  is  well  to  realize  that,  for  a 
specified  aspect  ratio  *  as  t^ie  l°a<*  is  increased  values 

of  A  corresponding  to  flexural  bifurcations  of  all  orders  are 
reached  before  any  barreling  bifurcation  is  attained.  Even  if 
these  could  be  inhibited,  the  barreling  bifurcation  of  highest 
(theoretically  infinite)  order  is  reached  first,  corresponding 
to  wrinkling  of  the  free  surfaces. 
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9. 


Appendix 
The  quantities 


and  £  introduced  in  §6  are  defined  by 


|  ra2(X2-l){(X-l)F3(X)  +  (X+1)F3(-X)}  , 

=•  4(16 [X2  (X2-l)m2  +  X2(X4+m4)  +  (X2+m2)2]  (9.1) 

-2Xm2(X2-l) [ (X+l) 3  +  (X-l)3] 

+  2 [  (X-l) 4  +  (X+l)4] (2X2+l)m2  -(X2-l)2(X2+l)m2}  . 


The  quantities  Cu(j,k)  and  ?5(j,k)  introduced  in  §6  are 
defined  by 

(j  »k) 

Cu(j,k)  =  fXTlc —  sinh{(jX+k)n> 

4  k  (9.2) 

?c(j  »k) 

C5(j,k)  =  JX+Tc -  sinh{  (j X+k) n) 

where  j  and  k  are  integers  and 


C4(3,l) 

CU(1.3) 

?u(2,2) 

Cu(2,-2) 

Cu(2,0) 

?u(0,2) 

?U(1,D 

?U(1,-D 


-X (X-l) 2mF  ,  ?u(3,-l)  =  -vX(X+l)2mF1  , 

-(X-l)2mF2  ,  ?u(l,-3)  =  -v(X+l)2raF2  , 


\  (X+l) (X-l)2m2F3(X)  , 

=  \  (X-l)  (X  +  l)2m2F3(-X)  , 

^v{16X3F1+  m2[(X+l)3F3(X)  +  (X-l) 3F3 (-X) ] } 
£v(16X2F2+  m2[(X+l)3F3(X)  +  (X-l) 3F3 (-X) ] } 
-v(X+l)m(4X2F3(X)  +  (X+l) (XF1+F2)}  , 

*  -(X-l)m{4X2F3(-X)  +  (X-1)(XF1+F2)}  , 


(9.3) 
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C  (4,0)  *  64AU(2A2-1)  ,  ?5C0,4)  -  64A2mU  , 

^(3,1)  »  -32A2 (A+l)m{ (A+l) (3A2-1)  +  2A(A-1)}  , 

?5(3,-l)  «  -32vA2(A-l)m{(A-l)(3A2-l)  -  2A(A+1)}  , 
f  (1,3)  *  -64A(A+l)m3(A2+2A-l)  , 

£^(1,-3)  *  -64vA ( A -l)m3 (A2- 2 A - 1)  , 

?5(2,2)  *  2m2{ (A-1)U[4- (A+l) (3A-1)] 

8A(A+l)3(A2+3A-2)  ♦  32A2  (3A2-1) }  , 

?5(2,-2)  *  2m2 { (A+1)^[4-(A-1) (3A+1) ] 

♦  8A(A-l)3(A2-3A-2)  +  32A2(3A2-1)}  ,  (9.4) 

^.(2,0)  -  4v(A2-l)(m2(A2-l) (5A2+3)  -  32 A2 (A2+m2) }  , 

?5(0,2)  -  4v(A2-l)m2{ (A2-l)  (3A2+5)  -  64A2}  , 

^  (1,1)  *  -32v (A-l)m{ (A2+m2) (A-l)  (A2-A+2) 

-  2A(A  +  1)  [m2+A(A2+-3A  +  l)]  }  , 

5  (1,-1)  »  -32(A+l)m{(A2+m2) (A+l) (A2+A+2) 

+  2  A ( A - 1 ) [m2-A(A2-3A+l)] }  . 
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Figure  Captions 

Fig.  1  Relation  between  n  and  critical  values  of  X  . 

Curve  I:  flexure;  Curve  II:  barreling;  Curve  III: 
asymptotic  expression,  correct  to  Ofn^)  . 

Fig.  2  Relation  between  G*  and  n  for  flexure,  with  X3  *  1  . 

Curve  I:  exact  expression;  Curve  II:  asymptotic  expression 
correct  to  0(n2)  . 

Fig.  3  Relation  between  G*  and  n  for  flexure,  with  X^  =  1  . 

Fig.  4  Relation  between  G*  and  n  for  barreling,  with  X^  * 
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of  an  energy  criterion,  without  restriction  on  the  thickness  of  the  plate. 

The  critical  state  corresponding  to  the  lowest  order  flexural  mode  is 
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found  to  be  stable  (unstable)  if  the  aspect  ratio  (thickness/length)  is 
less  (greater)  than  about  0.2.  Agreement  with  the  classical  Euler  theory 
is  established  in  the  limiting  case  of  small  aspect  ratio. 


